Insensitivity of sub-Kelvin electron-phonon coupling to substrate properties 
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We have examined the role of the substrate on electron-phonon coupling in normal metal films 
of Mn-doped Al at temperatures below 1 K. Normal metal-insulator-superconductor junctions were 
used to measure the electron temperature in the films as a function of Joule heating power and 
phonon temperature. Theory suggests that the distribution of phonons available for interaction 
with electrons in metal films may depend on the acoustic properties of the substrate, namely, that 
the electron-phonon coupling constant E would be larger on the substrate with smaller sound speed. 
In contrast, our results indicate that within experimental error (typically ±f0 %), E is unchanged 
among the two acoustically distinct substrates used in our investigation. 



At temperatures below 1 K, the interaction between 
electrons and phonons in metals weakens to the point 
that the two systems can be out of thermal equilibrium 
with one another. Such weakened coupling may be unde- 
sirable e.g., in dc superconducting quantum interference 
devices (SQUIDs) where the noise may be limited by hot- 
electron effects in the normal- metal shunts [l[ . This effect 
is also the origin of the problem of low-temperature ther- 
mometry wherein the electrons of the thermometer are 
hotter than the surrounding material. On the other hand, 
such decoupling may be desirable to increase the per- 
formance of electron-tunneling microrefrigerators based 
on normal metal-insulator-superconductor (NIS) junc- 
tions. Electron-phonon decoupling can also be exploited 
to make cryogenic bolometers where dynamic range and 
sensitivity are strong functions of the decoupling level Q ■ 

Electron-phonon (ep) coupling occurs because passing 
phonons distort the local lattice structure, and conduc- 
tion electrons respond to the resulting band distortions. 
A useful approximation for this interaction is the scalar 
deformation potential Q, which relates the change in 
Fermi level Sef to the local volume change or dilation, 
AV/V, of a unit cell, 
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For a longitudinal phonon, the dilation magnitude 
|AV|/V = y/hq/2pflc, where Ti is the reduced Planck's 
constant, q the phonon wavevector, p the mass density, 
n the crystal volume, and c the sound speed. A more re- 
fined theory of ep coupling, accounting for interaction 
between electrons and transverse phonons, was devel- 
oped by Reizer Later, Sergeev and Mitin [B| ex- 
tended Rcizcr's work by accounting for quantum inter- 
ference that occurs when static impurities are present in 
the host metal. Regardless of the framework employed, 
the two important quantities derived in these theories 
are the electron relaxation rate r^ 1 , and the heat flow, 
or thermopower, between the electron and phonon sub- 
systems P ep - 

Calculation of P ep involves integration over the elec- 
tron states k and phonon states q of an electron-phonon 



collision term M. t 
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where T e (T p ) is the electron (phonon) temperature. 
Converting to integration over energies w(q) and e(k), 
we have 
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where T>(u>) is the phonon density of states. Since A4 ep 
contains factors of the form / e (l — / e ), where / e (e, T e ) is 
the Fermi-Dirac distribution function, only the electron 
density of states at the Fermi level v(ef) is relevant, and 
this can be brought out of the e integral. Because of the 
factors v(ef) and T)(uj), as well as deformation potential- 
like terms in M ep , one finds that P ep is principally de- 
pendent on the Fermi velocity vf and the sound speed c 
of the metal. 

For bulk metals the ep theory outlined above is 
straightforward enough. However, many experimental 
situations that rely on ep coupling involve the use of thin 
metal films on dielectric substrates, and the correct form 
to use for T>(w) in Eq. <|3j> is not always clear. At low 
temperatures it is possible to be in a regime such that 
the thickness d of a film is much less than the dominant 
thermal phonon wavelength A p h ~ hc/AhsT, where fee 
is Boltzmann's constant, and the 4 in the denominator is 
based on an empirical conversion factor found in Ref. @. 
For a free-standing film, with A p h 3> d, the phonon dis- 
tribution becomes two-dimensional (2D) and the phonon 
density of states changes to reflect this new dimension- 
ality, as well as the presence of flcxural Raylcigh-Lamb 
modes 0]. On the other hand, if the film is strongly 
adhered to a rigid substrate, then the proper system to 
consider is that of the composite film/substrate, which 
should exhibit a three-dimensional (3D) phonon spec- 
trum. Further, for A p h ~> d it may be better justified to 
characterize that 3D phonon distribution with substrate 
parameters (e.g., p S ub and c su b), rather than those of the 
metal film. We report in this Letter the dependence of 
ep coupling on substrate, by comparing the measured ep 
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thermopower in normal metal thin films of Al doped with 
4200 ppm Mn fgj on two substrates (Si and Si02) with 
very different sound speeds (cf. Table HJ) . 

TABLE I. Longitudinal and transverse phonon speeds in [m/s] 
for materials used in this work. Also shown is the expected 
trend for E, based on electronic parameters of the Al:Mn film 
and sound speeds for each substrate. 



substrate 


longitudinal 


transverse 


expected £ 


Al 


6380 


3070 


Sai 


Si0 2 


5840 


3700 


~ Sai 


Si 


9420 


5860 





In order to quantify our results, we assumed the fol- 
lowing form for the ep thermopower: 

P ep = Eft (T e " - T p ") , (4) 

where fl is the normal-metal volume, n = 4 to 6, and 
S is the material-dependent ep coupling constant of or- 
der 1 nW//im 3 K" . This expression for P ep follows from 
completing the integration in Eq. ([3]). Depending on the 
particular regime (e.g., amount of disorder, temperature) 
theory predicts different expressions for S and different 
values for n [H, HI . The different regimes can be parame- 
terized by the product of the dominant thermal phonon 
wavevector = 4k bT '/he and the elastic electron mean 
free path £, as depicted in Table |TTJ 



TABLE II. Electron-phonon coupling regimes in bulk normal 
metals. 



qrl —¥ oo 


pure limit 


n — 5 


qre > 1 


clean limit 


4 < n < 5 


q T i < 1 


impure/dirty limit 


4 < n < 6 



The variation in n for qx£ ~ 1 is due to quantum in- 
terference between "pure" electron-phonon and electron- 
boundary /impurity scattering. In the clean limit, n — > 4 
if ep interaction is dominated by transverse phonons and 
scatterers have a static character — e.g., massive impu- 
rities or rigid boundaries. The same behavior is found 
in the dirty limit as the static component of impurities 
increases. In other words, an increase in the fraction S of 
static scatterers enhances ep coupling (n — > 4), while an 
increase in the overall concentration of scatterers weak- 
ens ep coupling [n — > 5 or 6). The theory developed in 
Ref . [5j for the dirty limit yields an expression for P ep that 
is actually the sum of n = 4 and n = 6 terms, 

P ep = £ 4 fi (T 4 - T p 4 ) + £ 6 ft (T e 6 - T p 6 ) . (5) 

As a result, when we fit the data to Eq. ((U we will obtain 
a single value for n that is between 4 and 6, depending 
on the relative contribution of the two terms in Eq. 
A similar expression exists for the clean limit, but with 



n = 5 instead of 6 and a different formula for £4. The 
reason for fitting the data to Eq. ([4]), and not Eq. ([5]) is 
that the theoretical expressions for £„ apply only to their 
respective limits, and as we discuss later, our data cover 
the crossover between the impure and clean regimes. The 
full expressions for £„ are complicated, but they all de- 
pend on the sound speed c as 



Therefore, if the substrate affects the phonon density of 
states 2?(w) in the metal film, then the measured E ra 
should differ by a factor in the range of 4-10 on the two 
substrates we have chosen (cf. Table Q] and Eqs. ([5]) and 
©)• 

In order to measure S, we fabricated thin- film resis- 
tors (thickness d = 35 nm) with integrated NIS (nor- 
mal metal-insulator-superconductor) thermometers (see 
inset of Fig. [1}, using standard photolithographic pro- 
cessing For our Al:Mn films and the temperature 
range covered, the dominant transverse phonon wave- 
length A p h = 60-370 nm, which is always greater than 
d. NIS junctions provide reliable thermometry as long 
as the product of the normal-state tunnel resistance Rt 
and junction area is greater than about 10 kf2-/im 2 . For 
these high-resistance junctions, normal-metal self-cooling 
is negligible, and the local electron temperature near the 
thermometer junction is identical to that in the rest of 
the normal-metal film. Heater current was injected into 
the resistor via superconducting electrodes. The result- 
ing NS contacts provide very good electrical conductance, 
but very poor thermal conductance, since they are always 
biased within the superconducting gap Aai- 

Since we are ultimately interested in determining the 
parameters S and n, and since the temperature depen- 
dence of ep coupling is roughly T 5 , accurate determina- 
tion of temperature is critical. The current-voltage (/V) 
curves of our NIS thermometers were found to fit rather 
well (cf. Fig[l} to the standard isothermal theory, based 
on the BCS density of states [lo| . 

(7) 

where e is the elementary charge and f(E) is the Fcrmi- 
Dirac distribution. Using the measured energy gap of Al 
(Aai ~ 185 /xeV) and the measured tunneling resistance 
Rt ~ 6 kfi (per junction) as inputs, and electron temper- 
ature T e as a fitting parameter, we fit IV data to Eq. (|7|) 
by use of a nonlinear least-squares algorithm. A loga- 
rithmic weighting was applied to the current (I) data, 
since the region of the IV curve below the gap has the 
greatest temperature response. This procedure allowed 
us to avoid some of the calibration difficulties associated 
with resistance thermometry. 

The IV curve fit at zero power determines the bath 
temperature Tbath, which is nearly equal to the cryostat 
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FIG. 1. NIS thermometer IV data (symbols) and correspond- 
ing fits (lines) for Tbath ~ 0.1 K. Black circles represent the 
unheated (P = 0) data, while red crosses indicate the data 
for an applied heater power P = 740 pW. The temperatures 
of the two fits are T — 118 mK (P = 0) and T = 275 mK 
(P = 740 pW). Inset: Schematic of NIS thermometer/heater 
device and measurement arrangement. Red resistor strip is 
Al:Mn. Grey current injection electrodes are superconducting 
Al. Resistor length is 500 /jm. 



temperature, as expected when stray loads are negligi- 
ble. By sweeping over all bath temperatures and applied 
powers P, and fitting the respective IV curves, we ob- 
tained a data set T e (P, Tbath)- The phonon temperature 
T p in the film was calculated from a model for power flow 
out of the phonon system that includes the predicted 



Kapitza resistance between the film and substrate 11 1 
and the predicted thermal conductivity of the substrate 
itself. Thermal resistance between the substrate and cop- 
per mounting box was not included in this model, but it 
was estimated to be much smaller than that between film 
and substrate. Ultimately, the calculated differences be- 
tween T p and Tbath are very small, so our final results 
arc insensitive to the details of this model. The final 
data set T e (P, T p ), along with our measurement of the 
heater strip volume, were used to determine E and n by 
fitting to Eq. (gj. 

The primary result of this Letter is presented in 
Fig. [2ja). Aside from a slight deviation among the two 
substrates at higher temperatures, it is readily appar- 
ent that changing substrates has a negligible effect on 
the ep coupling constant. The overall magnitude of E 
is consistent with previous measurements on Al:Mn thin 
films [l2|. The variation in E with temperature likely 
indicates a crossover from the impure limit at low T to 
the clean limit at high T. For example, assuming qr is 
a function only of the transverse sound speed in Al, we 
find 0.2 < qxl < 1.4 over the temperature range 0.1 K 
to 0.6 K for films on Si. Similarly, 0.4 < qx£ < 2.6 for 
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FIG. 2. (a) Measured £ and (b) n vs. electron temperature T e 
for Al:Mn resistors on Si (red, open symbols) and Si02 (black, 
filled symbols) substrates. The dashed and solid lines for each 
color represent different samples from the same wafer. Error 
bars indicate ± one standard deviation in the uncertainty of 
the fit. 



films on SiC>2. For the temperature range explored here, 
transverse phonons dominate the interaction with elec- 
trons. The characteristic temperature T* for transverse- 
longitudinal crossover can be estimated by comparing the 
electron relaxation rates of each polarization [j| in the 
regime qx£ ~ 1, 



T* = 



67T TlCl 

7C(3) k B £ 



(8) 



where £(3) ~ 1.202 and c; (c*) is the longitudinal (trans- 
verse) sound speed. For Al:Mn films on Si, T* ps 70 K. 
Because we are primarily concerned with the sub-kelvin 
regime, we can safely assume that transverse phonons are 
much more important for ep coupling. 

The clean-impure crossover is made more clear when 
considering the variation of the ep exponent n with T e , 
as shown in Fig. HJb). Here, it would seem that the sub- 
strate does exhibit an influence, but we hypothesize that 
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the difference in n between the two substrates can be at- 
tributed to a difference in the fraction of static scatterers 
where C is the mean free path due to static 
scattering. If we assume that scattering at the upper and 
lower film surfaces is predominantly static (l3j . the film 
thickness d sets the scale for C, and 6 ~ i/d. For our 
films on SiC>2 substrates, ^sio 2 ~ 25 nm, which is nearly 
twice that of our films on Si (£si ~ 14 nm) [sj]. Since d 
is the same for both films, <5sio 2 ~ 2<5si, and we expect 
that Al:Mn films on Si02 should indeed exhibit a smaller 
exponent than the films on Si, at a given temperature. 

The insensitivity of the ep coupling to the substrates 
studied here presents a puzzling problem. If the phonons 
of the substrate-film system are pictured as a gas that 
is able to move freely between the film and substrate, 
then phonons from the substrate will dominate the in- 
stantaneous phonon population in the film, since each 
body contributes a number of phonons to the gas that is 
proportional to its volume 0]. In this picture, substrate 
properties are expected to strongly influence ep coupling. 

Since £ may increase or decrease with £ Q, one may 
argue that the dependence of £ on c (i.e., the substrate) 
is hidden by virtue of the films exhibiting different £s 
on the two substrates. However, for T < 0.3 K, the 
n = 6 term in Eq. ([5]) is dominant for films on either 
substrate. For this case, S 6 ~ £/c 5 , and we estimate 
E6,Si0 2 /^6,Si ~ 18- In other words, at low temperatures 
the difference in £ should only enhance the difference in £ 
due to c. Since the observed Ss for the two substrates are 
practically equal for all temperatures considered here, it 
is reasonable to assume that the different mean free paths 
do not impact our overall conclusion that ep coupling in 
metal films is not strongly dependent on substrate. 

Decoupling of the phonon systems in the film and sub- 
strate could reduce the dependence of the ep interac- 
tion on the substrate. For example, phonon trapping 
due to acoustic mismatch could lead to modification of 
the phonon gas dimensionality in the film, despite the 
supported geometry used in this work. If the sound 
speed in the metal film is less than that in the substrate, 
phonons in the film experience total internal reflection 
when their angle of incidence at the film-substrate bound- 
ary is greater than a critical angle 9 C = sin _1 (cfii m /c su b) 
from the film-substrate interface normal This ef- 

fect is the acoustic analogue to Snell's law in optics. 
Only phonons that approach the substrate interface from 
within the cone are able to travel across the interface. 
Phonons traveling on paths outside the critical cone 
are trapped within the film and resemble a 2D phonon 
gas [ill. For Al:Mn on Si (SiC^), the angle of the critical 
cone for the dominant transverse phonons is 32° (56°). In 
other words, in films with a narrower critical cone (e.g., 
Al:Mn on Si), phonon trapping is enhanced and the film 
phonon spectrum may become 2D. However, for real in- 
terfaces other phonon scattering mechanisms exist, which 
tend to make the phonon gas 3D 15 1. 



Although there are many factors (e.g., temperature, 
disorder) that determine the observed ep exponent n, we 
can still use our data for n to help elucidate whether 
phonon dimensionality plays a role in our null result. 
First, the theory of Sergeev et al. @ — appropriate for 
bulk samples exhibiting a 3D phonon spectrum — pre- 
dicts that our Al:Mn samples should exhibit an n that 
varies from 6 to 4, as temperature is increased. This pre- 
diction is supported by our ep coupling measurements 
and those on similar Mn-doped Al films 12j. Second, it 
was indicated both theoretically and experimentally 
that ep coupling in the presence of a 2D phonon gas is 
associated with a value of n < 4.5. Finally, it may be 
argued that A p h is close enough to d that a crossover 
from 3D to 2D phonons occurs as temperature is low- 
ered. However, the observed trend for Al:Mn films on 
both substrates is for n to increase as T is decreased. 
If the film phonons become more 2D as T is lowered, 
we would expect n to decrease or remain fairly constant. 
These points strongly suggest that there is little meaning- 
ful difference in phonon dimensionality between Al:Mn 
films on the two substrates and that phonons in the films 
are closer to 3D than 2D. 

It is interesting that we obtain qualitative agreement 
with the theory of Sergeev et al. [E| by assuming the films 
behave as bulk samples (Eg = 4-14 nW//im 3 K 6 ). How- 
ever, since A p h is always greater than the film thickness 
d, it is not possible to construct phonon modes in the 
direction perpendicular to the film, and therefore yield 
a 3D phonon gas, without also considering the substrate 
as part of the overall problem. Our observation that 
the phonons are 3D, and yet ep coupling is insensitive to 
substrate, cannot be accounted for within this theoretical 
framework. 

In summary, we have demonstrated that ep coupling 
in thin normal-metal films is largely independent of the 
acoustic properties of the supporting substrate. The re- 
sults presented have important implications for efforts 
to engineer low-temperature micro- or nanostructures, 
which depend on ep coupling for operation. Further- 
more, the ubiquity of film-substrate geometries points to 
the need for models of ep coupling that consider such 
geometries explicitly. 

This work was supported by a National Research 
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copyright in the United States. 



[1] F. C. Wellstood, Excess noise in the dc SQUID; 4.2K to 
20mK, Ph.D. thesis, University of California, Berkeley 
(1988). 

[2] J. Wei, D. Olaya, B. S. Karasik, S. V. Pereversev, A. V. 
Sergeev, and M. E. Gershenson, Nature Nanotech. 3, 
496 (2008). 



5 



[3] C. Kittel, Introduction to Solid State Physics, 7th ed. 

(John Wiley & Sons, 1996). 
[4] M. Reizer, Phys. Rev. B 40, 5411 (1989). 
[5] A. V. Sergeev and V. Mitin, 

|Phys. Rev. B 61, 6041 (2000)| 
[6] T. Klitsner andR. O. Pohl, Phys. Rev. B 36, 6551 (1987). 
[7] See J. T. Karvonen and I. J. Maasilta, 

and references 



Phys. Rev. Lett. 99, 145503 (2007) 



[8] 



therein. 

G. C. O'Neil, D. R. Schmidt, N. A. Tomlin, 
Ullom, J. Appl. Phys. 107, 093903 (2010). 
[9] See Supplemental Material 

URL_will_be_inserted_by .publisher for 



and J. N. 

at 

additional 



details related to device fabrication and experimental 



setup. 

[10] M. Tinkham, Introduction to Superconductivity, 2nd ed. 
(Dover, 2004). 

[11] E. T. Swartz and R. O. Pohl, Rev. Mod. Phys. 61, 605 
(1989). 

[12] L. J. Taskinen and I. J. Maasilta, 
|Appl. Phys. Lett. 89, 143511 (2006)| 

[13] Static scattering at film boundaries may result from lat- 
tice distortions due to strain or reconstruction at the 
film surface and thermal expansion mismatch at the film- 
substrate interface. 

[14] S. B. Kaplan, J. Low Temp. Phys. 37, 343 (1979). 

[15] J. C. Nabity and M. N. Wybourne, Phys. Rev. B 44, 
8990 (1991). 



